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(N : Resumen 

En este arti'culo se presentan las ecuaciones de Euler y de Navier-Stokes, las mas 
basicas de las ecuaciones de los fluidos incompresibles, asf como unos modelos sim- 
plificados de dichos problemas como pueden ser la ecuacion quasigeostrofica u otros 
escalares activos. Las tecnicas matematicas para obtener que estas ecuaciones estan 
£f) . bien puestas en el sentido de Hadamard y otras propiedades cualitativas (Principios 

del Maximo, formation de singularidades...) se ilustran en el caso sensiblemente mas 
sencillo de la ecuacion de Burgers con disipacion no-local. Se adjunta una section de- 
dicada a los metodos numericos usados para aproximar soluciones a estos problemas. 
Este trabajo tiene su origen en una serie de clases que impartimos en la escuela JAE- 



o 

Intro del CSIC durante el verano del curso 2010-2011 y por lo tanto se centra en la 



X 



popularization y divulgation de la fisica involucrada y en la explication detallada de 
las ideas mas abstractas en los argumentos puramente matematicos. 
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1. Introduction y motivacion 

El estudio de las ecuaciones de los fluidos incompresibles tiene cada vez un mayor 
interes, tanto desde el punto de vista mas teorico (integrales singulares...) como desde el 
enfoque mas aplicado (simulaciones numericas...). 

Las ecuaciones que aparecen modelizando problemas de mecanica de fluidos son varia- 
das, pero las mas importantes son las de Euler y Navier-Stokes. De hecho por demostrar 
(o refutar) la existencia global de solution clasica para Navier-Stokes el Instituto Clay 
otorga un premio de un millon de dolares. 

Consideramos como dominio espacial el piano o el espacio enteros, i.e. M rf con d = 2, 3 
y exigimos que la velocidad del fluido sea cero en el infinito. Sea ademas / el campo de 
velocidades initial. Entonces las ecuaciones de Euler (1707-1783) para la velocidad u de 
un fluido incompresible son 

Masa Acclcracion Fuerzas internas Fuerzas externas 

p (dfU + (u ■ V)u)= — Vp + F (Conservation del momento), 

V ■ u = 0, (Conservacion de la masa). 

con V = (d xi , ...,d Xd ) y u = (ui, ...,Ud). Este sistema de ecuaciones es la segunda ley de 
Newton en el caso de un continuo de parti'culas. 

Las ecuaciones de Navier (1785-1836) y Stokes (1819-1903) para un fluido incompresible 
con densidad p son 

Masa Aceleracion Fuerzas internas Fuerzas externas 

/ ~ / p~ s {dtu + (u ■ V)it)= — + v Au + F (Conservation del momento), 

V • u = 0, (Conservacion de la masa), 

(2) 

donde A = 

En las ecuaciones de Navier-Stokes se ha anadido el rozamiento entre particulas del 
fluido modelizandolo con un laplaciano. 

Observamos que es un sistema de evolution 'extrano' en el sentido de que la derivada 
temporal de la presion p no aparece. Eso nos indica que la presion se puede obtener de la 
velocidad u. Para ver esto basta tomar la divergencia de la ecuacion de conservacion del 
momento, 

V • ((n • X7)u) = -Ap. 

Ahora podemos utilizar la funcion de Green para el Laplaciano (que para M. d es conocida) 
y obtener p = G(u). Tomamos el gradiente en esta expresion y obtenemos un sistema 
de ecuaciones no-locales cerrado para u (donde ahora nos restringimos a velocidades in- 
compresibles). Ademas hemos obtenido que la presion actua como un multiplicador de 
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Lagrange para la restriction V • u = 0. Los resultados de los que se dispone cierran la 
teoria en d = 2. En el piano se conoce la existencia global (en tiempo) de soluciones 
clasicas. Para el problema completo d = 3 se tiene un teorema de existencia local (en 
tiempo) de solution clasica (ver [19]). Esta diferencia es porque las ecuaciones en los casos 
bidimensional y tridimensional son sensiblemente distintas. Para verlo necesitamos defi- 
nir una cantidad que nos de information sobre cuanto 'gira' el fluido. Esta cantidad que 
denotaremos u se llama vorticidad y es el rotational de la velocidad, i.e. 



con u el campo de velocidades del fluido. En concreto se tiene que las ecuaciones para la 
vorticidad son distintas: para un fluido no viscoso e incompresible 1 , es decir, que sigue las 
ecuaciones de Euler, en el caso d = 2 se tiene 



que es una ecuacion de transporte, mientras que en el caso d = 3 se tienen las siguientes 
ecuaciones 



En el caso de consider ar un fluido viscoso se ha de anadir un termino Aw. 

La diferencia entre ambos casos aparece tambien en las longitudes de onda que estan 
mtimamente relacionadas con el fenomeno de la 'turbulencia' (ver [13]). La turbulencia 
tiene como efecto principal facilitar que dos fluidos se mezclen, por lo tanto, llegados a 
este punto podemos 'experimentar' un teorema. Para este pequeno juego necesitamos dos 
vasos pequenos Uenos hasta arriba uno de ellos de agua y el otro de vino. La cuestion es: 
icomo conseguimos cambiar los Kquidos de vaso sin usar un tercer recipiente y sin que se 
mezclen? Para responder a esta pregunta hemos de conocer como es la turbulencia en tres 
dimensiones y que diferencia hay con dos dimensiones. Asf, si conseguimos una manera 
de reducir el problema tridimensional a uno bidimensional hemos acabado, porque en dos 
dimensiones 'no hay turbulencia' y esta es la culpable de que los hquidos se mezclen. Para 
conseguir esta reduction en las dimensiones lo que hacemos es tapar el vaso de agua con 
un carne y colocarlo con cuidado encima del vaso de vino. Si lo hemos hecho bien no se 
ha salido ni una gota. Ahora abrimos una rendija minuscula entre los vasos y el carne. El 
agua es mas densa, por lo tanto comenzara a bajar mientras que el vino subira...jy todo 
esto sin mezclarse! (ver Figura 1) 

Tras esta excursion por las ciencias experimentales volvamos a las matematicas. La 
vorticidad u es una cantidad que aparece en el conocido criterio de existencia global de 
Beale-Kato-Majda (ver [4]). 

Teorema 1 (Beale-Kato-Majda). Sean T* y M dos constantes tales que 



entonces una solucion clasica de las ecuaciones de Euler (1) existe al menos hasta tiempo 
T* . Ademds, si T max es el tiempo mdximo de existencia (es decir, aparece una singulari- 
dad), entonces 



rot u = co 



dtu + u ■ Vw = 



(3) 



dtu; + (u • V)cj = (co • V)«. 



(4) 





1 Estos fluidos se conocen como fluidos ideales. 
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Figura 1: Experimento. 



Este teorema nos dice que si controlamos 'lo que gira' el fluido entonces no hay sin- 
gularidades, por lo tanto conocer como se comporta la vorticidad es crucial para intentar 
entender que hace el fluido. En efecto, tambien es interesante porque si la conocemos 
podemos recuperar la velocidad gracias a la formula de Biot-Savart (ver [3]). En el caso 
bidimensional la ley de Biot-Savart es 



Ahora bien, con lo que acabamos de mencionar podemos ver que el problema de las 
ecuaciones de Euler (1) en el caso bidimensional podemos formularlo de manera equivalente 
como un 'escalar activo', es decir, un escalar que es transportado por el fluido de manera 
que ademas podemos recuperar la velocidad del fluido si conocemos el escalar. En efecto: 
recordemos que temamos la ecuacion para la vorticidad (que en dos dimensiones es un 
escalar) (3), esto unido a u = K{uj) cierra el problema para la vorticidad. Ademas si 
asumimos que nuestro dato inicial esta acotado, usando el Teorema 1, tenemos que la 
solution clasica del sistema (1) existe globalmente. 

El estudio matematico de los escalares activos tiene una gran relevancia en cuanto 
que son sistemas sencillos que conservan el caracter no local de un fluido incompresible. 
Ademas de que hay multitud de aplicaciones donde el problema fisico se puede modelizar 
con un escalar activo. Por ejemplo tenemos el caso de la ecuacion Quasigeostrofica (ver 
[5] j [8], [16], [7] y las referencias allf expuestas) o de la Ley de Darcy (ver [10], [11], [12], 
[9] j [1], [20]). La ecuacion quasigeostrofica modeliza la evolution de la temperatura de 
grandes masas de aire en grandes escalas espaciales y es muy estudiada como modelo de 
la 'frontogenesis' (la formation de frentes de aire a distinta temperatura). Este problema 
es de interes en meteorologfa, porque ya se sabe que 




(5) 



con 




La falta de acierto de quienes predicen el tiempo se ha hecho ya proverbial, 
y sin embargo no hay ningiin meteorologo competente que no opine que los 
procesos atmosfericos estan causalmente deter minados. 

Max Planck (extrafdo de [21]) 

La ecuacion quasigeostrofica en dos dimensiones espaciales es 

d t 6 + u • V0 = 0, 
u = R ± u = (-R 2 9,R 1 6), 



(6) 



donde Ri es la transformada de Riesz i— esima (ver [22]) y 8 es la temperatura del aire. 

La Ley de Darcy modeliza un fluido incompresible que se mueve a bajas velocidades 
por un medio poroso. Asf si p(x, t) es la densidad del fluido se tiene el problema 

d t p + u ■ V p = 0, 
u = -(X7p + ge 2 p), 
V-it = 0. 

Este problema puede reducirse a un escalar activo tomando el rotacional dos veces en la 
ecuacion, obteniendo 

d t p + u-Vp = 0, . . 

u = (— A) _1 (rot rot g&iP)- 

Podemos modelizar la transferencia de calor interno del fluido en un medio poroso con 
la misma ecuacion (7) si a p le damos el sentido de una temperatura (ver [2]). Por lo 
tanto, tanto en el caso de la ecuacion (6) y de (7) puede interesarnos ahadir un termino 
de difusion del calor. Sin embargo estos terminos de difusion no tienen por que ser el 
tfpico laplaciano, puede ser necesario ahadir una 'potencia fraccionaria del laplaciano '. El 
laplaciano en el espacio de Fourier (con variables £) tiene una expresion sencilla, es un 
multiplicador: 

^Au = |£| 2 u. 



Podemos definir el operador A = \/—A de la siguiente manera 

&u=m, (8) 

y equivalentemente 

A^u = \£\ a u. (9) 

Notemos que tambien podemos escribir el resultado de aplicar el laplaciano fraccionario 
como la siguiente convolucion: 

A%(x) = /3(a,d)PV. / ^ ~ ^ dy (10) 
Jr>i \x - y\ d+a 

donde /3(a, d) es una constante de normalizacion. Los operadores (8) y (9) (o su version 
(10)) son lo que nosotros entendemos por 'potencias fraccionarias del laplaciano'. Asf las 
ecuaciones (6) y (7) con difusion no local son 

d t e + u- ve = --fA a e, 

u= R ± u = (-R 2 6,R 1 6), [ ' 



d t p + u-X7p = -^A a p, 
u = {R l R 2 ,-R 2 1 )p. 



2^ (12) 
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Para irnos aproximando a estos problemas podemos plantearnos otros mas sencillos o 
simplificados. Por ejemplo, consideremos un escalar activo 

J d t rj + u ■ Vtj = -^A a r], , . 

\u = T( V ), [i6) 

donde T es un operador integral singular que ademas nos garantiza que V-u = 0. Queremos 
simplificar este problema de manera que sea facilmente abordable, pero no tiene que estar 
tan simplificado que no nos de ninguna informacion. Lo primero que hacemos es reducir el 
mimero de variables espaciales a una, por lo tanto V = d x . Ademas podemos simplificar T 
tomandolo igual a la identidad, es decir, pierde su caracter integral singular. La ecuacion 
resultante de estas simplificaciones es la ecuacion de Burgers viscosa si 7 > y la ecuacion 
de Burgers no viscosa si 7 = (ver [17] y [14]): 

d t r] + nd x T] = --yA a n. (14) 

Como la condition V • u = en una dimension no tiene sentido fisico (en una dimension no 
hay choques) con esta ecuacion tenemos un modelo unidimensional de las ecuaciones (1) 
y (2). Asi nuestra ecuacion (14) es el modelo mas sencillo que nos da informacion tanto 
sobre las ecuaciones (1), (2) como (13). 

Este texto esta organizado de la siguiente manera: en la seccion 2 probaremos la exis- 
tencia local de solucion clasica para toda la familia de ecuaciones (14). En la seccion 3 
obtendremos una ley de conservation y unos principios del maximo que se tienen para las 
soluciones clasicas de la ecuacion (14). En la seccion 4 daremos un criterio de existencia 
de solucion clasica analogo al de Beale-Kato-Majda. En la seccion 5 veremos que hay blow 
up para 7 = en (14) y en la seccion 6 haremos simulaciones numericas de la solucion de 
(14). 



2. Existencia local de la solucion clasica 

El problema que estudiaremos en esta seccion y las siguientes es la ecuacion (14) (donde 
cambiamos la notation 77 por la mas corriente u): 

( d t u + ud x u = --/A a u, (t, x) G [0, T] xl 

\ u(Q,x) = f(x), {Lb) 

donde T > es el tiempo de existencia, 0<a<2y7>0. Sobre el dato inicial 
impondremos las condiciones necesarias cuando veamos el teorema de existencia local, 
pero por el momento podemos suponer que / £ C 2 . Nos restringiremos a las soluciones 
clasicas o soluciones fuertes. Este tipo de soluciones es tan suave como sea necesario para 
dar el sentido usual a las derivadas parciales que aparecen en la ecuacion. Asf para la 
ecuacion de Burgers queremos que u tenga al menos una derivada en i y dos en x. Es 
decir, queremos que u G . 

En esta seccion probaremos la existencia local de solucion clasica para la ecuacion (15) 
con 7 > y a G (0, 2] si el dato inicial esta en un cierto espacio de funciones. Para ello 
utilizaremos el metodo de la energfa de Leray (ver [19]). El argumento de Leray ya es un 
argumento clasico. Nosotros trataremos de aplicarlo con todo detalle en el ejemplo sencillo 
de la ecuacion (15). 

La idea del metodo es conseguir una sucesion de problemas regularizados para la ecua- 
cion (15). Para todos los problemas regularizados se demuestra la existencia utilizando el 
Teorema de Picard en un espacio de Banach adecuado. Se concluye el argumento obser- 
vando que la familia de soluciones regularizadas forma una sucesion de Cauchy, y por lo 
tanto convergente en algun espacio. El espacio de Banach que vamos a usar es H S (M), 
s > 3 porque utilizando la inmersion de Sobolev tenemos que entonces u G C 2 (R). 
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2.1. Estimaciones a priori 

Para utilizar el metodo de Leray hemos de conseguir unas cotas a priori para ciertas 
cantidades. Es decir, suponiendo que hay solution. 

Si pensamos en u como en la velocidad de un fluido incompresible nos interesa conocer 
que ocurre con la 'energta cinetica' del fluido, es decir, la norma I? de u. 

Proposicion 1 (Principio del maximo para || • Wi?). Sea u una solution cldsica de (15). 
Entonces se tiene que: 



a) Si 7 = 0, 

b) Sij> 0, 



W\\ L <t) = \\f\\ L 2, 

MU»(*)<ll/llr*. 



Demostracion. Observamos que 

Jm. Jm it 



1 d 



d x (u 3 )dx — 7 / A a uudx = — 7 / A a uudx. 



De aquf se concluye la parte a) del resultado. Para obtener la parte b) hemos de utilizar 
el Teorema de Plancherel: 

K a uudx = / \[\° illicit = I (A" -a) 2 ,lx. 



Jm. Jr ^ 



De esta ultima igualdad (que no es mas que la prueba de que el operador A a es simetrico) 
se concluye el resultado. □ 

En el caso de que la solution u no sea lo bastante regular como para ser solution 
clasica entonces la norma I? no se conserva. Esto en las ecuaciones de Euler es una serie 
de celebrados arti'culos (ver [18]). Las soluciones debiles de (1) en tres dimensiones o 
de (15) con 7 = no siempre tienen sentido ffsico. Lo que implica fisicamente es que 
un fluido perfecto que estuviese inicialmente en reposo puede comenzar a agitarse a lo 
loco sin haber mediado fuerza externa alguna, siempre y cuyo la velocidad u no tenga 
la suficiente regularidad como para ser solution clasica. Matematicamente esta 'paradoja' 
fisica se traduce en que no hay unicidad de soluciones debiles para el sistema de ecuaciones 
(1) en tres dimensiones. En el caso de (15) sin viscosidad (7 = 0) lo que ocurre es que se 
da un choque de curvas caracteristicas y despues no hay una manera unica de continuar la 
solution. Enlaza esto con las condiciones de Rankine-Hugoniot y las soluciones de entropia 
(ver [15]). Ademas, si consideramos el problema viscoso 

dtu 1 + v?d x u> = - 7 A°u 7 , 

pero tomamos 7 — > recuperamos la solution de entropia para la ecuacion de Burgers no 
viscosa como hmite, es decir, u° = hm 7 _>.o u 7 (ver Figura 2). 

Ya tenemos una cota a priori para la norma L 2 . Como mencionamos antes queremos 
que nuestra solution este en el espacio de Sobolev H 3 , por lo que solo falta estimar la 
norma I? de dtu. 
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Figura 2: Soluciones para distinas 7. 



Proposicion 2. Sea u una solution cldsica del problema (15) con un dato initial f G 
ff 3 (]R). Entonces se tiene la siguiente cota 



Demostracion. Se tiene que 



, } (]l 1 1 -x "Ilia - / u x uu x ufuux, 

por lo que hemos de derivar tres veces la ecuacion (15), obteniendo 

dld t u = -~/A a d%u - 3(d 2 x u) 2 + 4d x ud%u + u&^u. 
Introducimos esta expresion en nuestra expresion y logramos 

1 "'"d 3 u||| 2 = -7 / d*uA a dludx-3 [ {d 2 x ufdludx+ [ 4d x u(d^u) 2 dx 



\u\\ H s(t)< '"''j? (16) 

1 - C||/||#3t 



^-t-||9 3 m|| r2 = / diudidtudx, 



2dt 

+ / ud x ud x udx = h + I 2 + h + h- 



Usando el Teorema de Plancherel se obtiene que 

h = -7 / \t\ a dludludt = -7 / (k a ' 2 dlu) 2 dx = - 7 ||A Q/2 d x 3 U ||i 2 < 0. 

JR JR K 7 

La segunda integral se anula por las condiciones de borde impuestas 



I 



2 



[ d x [(d 2 x uf] = 0. 

JR. 



Las integrales que faltan son las mas singulares por tener el mayor niimero de derivadas, 
sin embargo solo nos hemos de preocupar de ^3, porque I4 es igual. En efecto 



h = \ [ ud x [(d 3 x u) 2 ] =-\( d x u{dlu) 2 . 

1 JR Z JR 



8 



Por lo tanto solo hemos de acotar la integral 1%. Se tiene que 

h < c||9 a! u||i, c||a^'u||^2 < \\u\\% 3 , 

y conchnmos que 

d 1 1 _Q 1 1 O II 1 1 Q 

^110*1*11X2 < C||«|| H 3. 

Usando la Proposicion 1 obtenemos 

d ii2 

Sin mas que integrar la EDO se concluye la primera parte del resultado. □ 
Comentario 1 Es mas, si / G H k (M.), k > 3 se tiene la siguiente cota 



\u\\ Hk (t)< Iff (17) 

1 - c\\f\\H"t 



2.2. Regular izacion del problema 

Consideremos una funcion p £ C£°, p(x) = p(\x\), p > y j R p = 1. Un ejemplo de tal 
/) puede ser la funcion 



p(x) 



f\e 1 -v 2 dy 

Ahora, para cualquier e positivo, consideramos las funciones 

Observamos que p e siguen teniendo las mismas propiedades que p. Dada una funcion 
g 6 L p , si ahora hacemos la convolucion g e = p e * g, obtenemos que <7 e G C°° y, gracias a 
la desigualdad de Young: 

111 

\\Pe * 9 \\L-r < \\Pe \\Li \\9 \\Lp, - + 1 = - + -, (19) 

r p q 

tenemos que 

\\Pe * g\\LP < \\9\\lp- 

Ademas lfm e _>o5e = 9- Pueden verse mas propiedades de estos micleos p e en [3]. 

El primer paso del metodo es regularizar el problema de manera que podamos demostar 
existencia de los problemas regularizados y tambien podamos usar las estimaciones a 
priori. Para ello utilizarmos la convolucion con p t . Asf consideramos la familia de problemas 

d t u e = -p e * (p e *u e )d x (p e *u e ) - p e * (jA a p e *u e ) , (t, x) G [0, T] x R,e > . . 
u e (0,x) = f(x). (ZU) 

Esta regularizacion particular se ha elegido pensando en obtener las mismas estima- 
ciones de la Proposicion 2.1. Por ser p e una funcion radial tenemos que 

fPe*9 = J J f(x)p e (x - y)g(y)dydx = J J f{x)Pe{y ~ x)g{y)dydx = j gp e * /, 

e integrando por partes tenemos unas estimaciones para u e como las de la Proposicion 2.1. 

Proposicion 3. Sea u e solution de (20) con dato initial f G H k (R), k > 3. Entonces se 
tiene la siguiente cota 



\u e \\ Hk (t)< , . . (21) 

1 - c\\f\\ H kt 
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2.3. Existencia para los problemas regularizados 

Para demostrar la existencia de los problemas regularizados (20) vamos a utilizar el 
teorema de Picard en el espacio H k (R). Hemos de ver entonces que para todo e se tiene 
que 

F e = -p e * (p € * U e ) d x (p e * U e ) - p e * {jA a p e * U e ) 

es localmente (es decir, si < A para cierto A) Lipschitz con respecto a la norma 

H k . Para empezar hemos de asegurarnos que no per demos derivadas al aplicar F e . Esto se 
consigue porque en lugar de derivar u t derivamos el micleo p e , de manera que u t no pierde 
derivadas. Veamos que F e es Lipschitz: 

Lema 1. Sea A > 0, entonces para todo e > F t es una funcion Lipschitz en {g : g G 
H k (R),\\g\\ Hk < A} con k > 3. 

Demostracion. Se tiene que 
\\F € (u) - F £ (v)\\ H k = 1 1 - p e * (pe * u) d x (p e *u) - p e * {-yA a p € * u) 

+ Pe* {p e * V) d x (p e * V) + p e * {jA a p € *v)\\ H k. 

Vamos a agrupar las partes difusivas de u con la de v: 

\\p e * (-jA a p t *u + jA a p e * v) \\ H k < tI |Pe *(v- u)\\ H k+ a < Li(7,a,e)||n - v\\ H k. 
Queda probar la estimacion para la parte convectiva: 

\\p e * (~p e * Ud x (p e *U) + p e * Vd x (p e * v))\\ H k 

< || - p e * Ud x (p e *u)±p e * Ud x (p e * v) + p t * vd x (p t * v)\\ H k 
< \\p e * ud x (p e *(v- u)) + (pe * (V - u))d x (p t * v)\\ H k 

< \\p e * Ud x (p t * (V - U)) \\ H 3 + \\(pe * (v - u))8 x (p t * v)\\ H k. 

Utilizamos ahora que, si s > 1/2, H S (R) es un algebra de Banach, es decir, que se cumple 

\\fg\\H* < ||/IMM|ijs. 

Gracias a esta propiedad de los espacios de Sobolev tenemos que 

\\pe * ud x (p t *(v- u))\\ H k + \\(p e * (v - u))d x (p e * v)\\ H k 

< \\p e * u\\ H k\\d x (pe *(v- u))\\ H k + \\d x (p e * v)\\ H k\\(p e * (v - u))\\ H k 

< L 2 (e, A)||u - v\\ H k. 

Para concluir hemos de elegir L = max{Li, L2}. □ 

Por lo tanto, dado un dato inicial / £ H k , si aplicamos el Teorema de Picard tenemos 
que existe una sucesion u t £ C' 1 ([0, T e ], H k (R)). Ademas, por las propiedades de los nucleos 
p e , se puede demostrar que T e = 00. 

Tenemos asf el siguiente result ado: 

Lema 2. Dado una dato inicial f G H k con k > 3 el problema (20) tiene una unica 
solucidn u e G C 1 ([0, 00), H k ). 
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2.4. Paso al h'mite e — > 

Veremos que u € G C 1 ([0, T], i^ fc ) forma una sucesion de Cauchy en el espacio C([0, T], if s 

< s < k donde T < T* = -njy, — (ver Proposicion 3). Para ello tenemos que ver que en 

I? es Cauchy y entonces utilizando las estimaciones de la Proposicion 3 concluiremos que 
es una sucesion de Cauchy en H s con s < k. El caso extremo de s = k lo trataremos al 
final. 

Lema 3. u e es una sucesion de Cauchy en C([0, T], H s ), con < s < k. 

Demostracion. Veremos primero que u e es Cauchy en C([0, T], I?). Sean e, 5 dos mimeros 
positivos. Entonces se tiene 

1 d f 

- — 1 1 u e -u s ||^ 2 = -7 / (p e * (A a p e *u € ) - p s * (A a p 5 * u s ))(u t - u 5 )dx 
l at Jro 

(p e * (p e * u e (d x (p £ * u e ))) - ps * (ps * u s {d x (ps * us))))(u e - us)dx 

= h + h- 

Veamos el caso de los operadores difusivos: 



h = -7 / (Pe * (A a p t * u e ) - ps* (A a ps * u s ))(u e - us)dx 
jr 

= -7 / (Pe * (A a p e * u e ) ±ps* (A a p 5 * u e ) - ps * (A Q p s * u s ))(u e - u 5 )dx 

JR 

= -7 / ((p e * (A a p e ) - ps* (A a ps)) *u t + ps* (A a p s * (u t - us))(u e - u s )dx. 

JR 

Integrando por partes el segundo sumando obtenemos que 

- 7 / ((Pe * (A>e) - Ps * (k a Ps)) *u € + p 5 * (A a p 5 * (u € - u s ))(u t - us)dx 

JR 



< -7 / (p e * (A a p t ) - ps* (A a p 5 )) * u e (u e - us) 

JR 

< j\\u t - us\\ L 2\\(p t * (A a p e ) — ps * (A a ps)) * u e \\ L 2, 
y por lo tanto, usando las propiedades de las suavizaciones (ver [3]) 
\\A a ((p e *(p t )- p 5 *(ps))*u e )\\ L 2 < \\(p t *(p t )- ps*(ps))*u t \\ H <* 

= \\(p e * (p e ) - ps* (ps)) *U e ±p e *U e ±p S * U e \\ H <* 

< \\p t * (p e * lt e ) - p t * U e \\ H a 
+ \\PS * (PS * U e ) ~ PS* Ue\\H<* 
+ \\Pe *U t - PS* U e \\H a 

< ce\\p e * + c5\\ps * 

+ \\pe * U e ± U e — PS * U € \\H a 

< Ce\\pe * Ue\\jjl+a + C(5||p5 * U e ||^l+« 

+(ce + c5)||n e ||^i+ Q 

< cmax{e,5}||« e ||^3 

Recordemos que, fijo T < T* , tenemos la cota uniforme ||u e ||#3 < C(T, WfWu^) Y P or 1° 
tanto tenemos que 

h < C(T,\\f\\ H 3,-f)\\ue - u s \\ L 2 max{e,<5}. 
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Para la integral I2 se hace igual. Hay que sumar y restar los terminos 

ps * (p e * u t d x (p € * u e )), ps * (ps * u € d x (p e * u e )), ps * (ps * u s d x (p e * u e ))... 

y aplicar las mismas desigualdades (Sobolev, Cauchy- Schwartz...). 

Una vez que hemos probado que u e es una sucesion de Cauchy en C([0,T],L 2 ), utili- 
zando la cota uniforme en e (Proposicion 3) 

\\u e \\ Hk <C(T, \\f\\ Hh ), 

podemos interpolar (ver [3]) entre el espacio H° = L 2 y H k , de manera que u e es de 
Cauchy para todo H s con s < k. Podemos asi asegurar la existencia de u G C([0, T], H s ) 
como el hmite de u e en la topologi'a usual de ese espacio. □ 

Ademas podemos mejorar el resultado para cubrir el caso extremo s = k. Asf hemos 
de ver que u = h'm e ^o u e esta en el espacio H k . 2 

Lema 4. La funcion u definida como el Umite anterior esta en C([0, T], H k ). 

Demostracion. Para probarlo hemos de razonar con analisis funcional. Como {u e } esta aco- 
tada en H k (Proposicion 3) tiene una subsucesion que converge debilmente en este espacio. 
Es decir 

< u e , z >-?■< v, z >, VzG H~ k . 

En principio no podemos asegurar que u = v. Para concluir la igualdad utilizamos las 
convergencias debiles, pues la convergencia fuerte implica la convergencia debil, en los 
espacios H s con s < k ademas del hecho de que H k C H s =^ H~ k D H~ s y obtenemos 
que si u 7^ v tenemos una contradiccion. □ 

2.5. Regular idad de u 

Gracias a los lemas anteriores podemos demostrar el siguiente resultado: 

Teorema 2 (Existencia local de solucion clasica). Sea f € H k , k > 3. Entonces existe una 
unica u G C([0, T], H k ) n C 1 ([0, T], H k ~ 2 ), definida como el Umite uniforme en compactos 
de R de u t , solucion clasica del problema (15) con 7>0|/0<a<2. 

Demostracion. Existencia: De los lemas anteriores se obtiene la existencia dew £ C([0, T],H k ) 
como hmite uniforme en compactos de u e . Para obtener la regularidad temporal necesa- 
ria para u observamos que dtu e tiende como distribution a la derivada debil de u (que 
denotamos por dtu). En efecto, para toda funcion test <p, se tiene que 

4>dtu e dtdx = — / dt4>u e dtdx — > — / / dt<fiudtdx = / / 4>dtudtdx. 
Jr Jo Jr Jo Jr Jo 

Ademas dtu t — > —~/A a u — ud x u en C([0, T], C(R)) (y por lo tanto tambien en sentido 
distribucional). Esto es una consecuencia de las propiedades de los nucleos suavizantes. 
Ahora estamos en condiciones de conseguir una mejor estimation para dtu: 

\\dtu\\ri\nT] rtu^ = max max I — jA a u — ud T u\ < max (c\ \u\ \ r 2 + I Inl |^ 2 ) 



< c\ 



\ u \\c([0,T],H 3 (R)) + \\ u \\c([0,T],H 3 (M.)) - C (\ l/l \n k i T , l)- 



Podemos elegir la notion de lmiite de cualquier espacio C([0, T], H s ) con s < k. 
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En particular si u G C([0,T], H k ) con k > 3 se concluye que -u es solution clasica de la 
ecuacion de Burgers. 

Para concluir que u 6 C 1 ([0, T], H k ~ 2 ) tenemos que obtener una cota para | \dtu\\ H k~2 : 
\\d t u\\ H k-2 < r y\\A a u\\ H k-2 + \\u\\ H k-2\\d x u\\ H k-2 < C(||/|| H fc,T,7), 

y entonces 

\\ u \\c 1 ([0,T],H k - 2 ) = maxJMItf*- 2 + PHI/^" 2 

- o^t"""^*" 2 + ™4T 7 " AQn " Hfe " 2 + r ?^ll n ll^ fe - 2 H 5:cn ll^ fe - 2 - cr (ll/llff fc » r >T0' 

donde en la ultima desigualdad hemos usado (17). 

Unicidad: Supongamos que hubiese dos soluciones, u y v en el espacio H k con k > 3 
del problema (15) con el mismo dato inicial. Entonces 



1 d .. 
2dt u 



|2 
L 2 



7 ||A Q / 2 (n-t;)||| 2 - / (u-^^tt-uS^di 



< -i / (u - v)d x (u 2 - v 2 )c!.r 

-1 
2 



< — - / (u — v)5 x ((u + «)(tt — v))dx 



< — - [ (u — v) 2 d x (u + v)dx — — [ (u — v)d x (u — v)(u + v)dx 

2 ii 2 Jr 

< c / (u — v) 2 d x (u + v)dx 

Jr 

< c(\\u\\ H k,\\v\\ H k)\\u-v\\ 2 L 2 

< c(T,\\f\\ Hk )\\u-v\\ 2 L 2, 

y usando Gronwall se obtiene la unicidad. □ 

3. Propiedades de las soluciones 

En esta section encontraremos algunas propiedades basicas que debera tener una so- 
lution de la ecuacion de Burgers. 

Como la ecuacion (15) modeliza un escalar que es transportado por un fluido donde la 
velocidad del fluido la hemos hecho proportional al escalar (14) es de esperar que la masa 
total, i.e. se conserva. En efecto, 

Proposicion 4 (Conservation de la masa). Sea u una solution clasica de la ecuacion (15) 
entonces se tiene que 

I u(t, x)dx = I f(x)dx. 
Jr Jr 

Demostracion. Integrando la ecuacion (15) en espacio 

— f u(t,x)dx= I dtu(t,x)dx= f — \-d x (u(t, x) 2 ) — jA a u(t, x)dx = I\ + 1^. 
dt Jr J Jr 2 

Por las condiciones de borde para u la integral I\ se anula. 
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Para ver que la integral I2 se anula razonamos utilizando la transformada de Fourier. 
Dada un funcion g{x) se cumple que 



g(x)dx = g(0). 

Como A a en el espacio de Fourier es el multiplicador de la definition (9) obtenemos que 

I 2 = - 7 / A a u (t,x)dx = -jA^u(0) = 0. 



□ 

Otra cantidad de la que es interesante ver la evolution es la norma L°° de la solution: 

Proposition 5 (Principio del maximo para 1 1 • | ). Sea u una solucion cldsica de (15) tal 
que max^gK u{t, x) > 0, rnin^gi u(t, x) < y a G (0,2], entonces se tienen los siguientes 
resultados 

a) Si 7 = 0, 

(t) = ll/l , 

b) Sij> 0, 

\\u\\ L °°(t) < H/Hloo. 

Demostracion. Observamos que la funcion M(t) = max^R u(t, x) = u(t, xt) es Lipschitz. 
En efecto: 

max \u(ti, x)\ = m&x(\u(ti, x)—u(t2, x)+u(t2, x)\) < max(|u(ii, x)— ufe, x)|)+max \u(t2,x)\, 

X X X X 

de donde 

I maxn(ti,x) — maxu(t2,x)| < max(|u(ii,x) — u(t2,x)\) = max(|<9 t n(s, x)\ \ti — ^ 2 1 ) 

XXX X 

< max max(|<9 4 ii(s, — ^2 1 - (22) 

se(t 2 ,ti) x 

Para concluir que la funcion es Lipschitz observamos que debemos restringirnos a un 
intervalo temporal [0, T] con T > fijo. Ahora nos aseguramos de que (ti,*2) C [0, T] y 
obtenemos 

I maxu(ti,x) — m&xu(t2,x)\ < max m&x(dtu(s, x))(ti — £2) = L(t\ — £2)- 

x x s€[0,T] x 

Usando el Teorema de Radamacher tenemos que M(t) es derivable en casi todo punto, 
M , = ^ M(t + hj)-M(t) = Km u(x t+h] ,t + h,)-u(x t ,t) 

hj-*0 hj hj-+0 hj 

u(x t+h t + hj) ±u(x t+hj ,t) - u(x t ,t) 
= bm J - J - = d x u(t, x t ) + o t u(xt, t) 

hj^Q hj 

= d t u(x t ,t). 
De manera que, si 7 = 0, 

M'(t) = d t u(x t ,t) = -u(x t ,t)d x u(x t ,t) = ~d x (u 2 (x t ,t)) = 
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donde en la ultima igualdad hemos usado que u(xt,t) es maximo. En el caso en el que 
7 > tenemos que dar signo al termino — r yA a u(t,xt). Usaremos la expresion como una 
convolution (10). Se tiene que 

- 7 A^.«)=eP.V./ ^»-ffi'U <0. 

jr f ~~ y\ 

En el caso en el que a = 2 el resultado se reduce a conocer que signo tiene la segunda 
derivada en un punto de maximo. 

Para el nrinimo se hace de manera analoga. □ 



4. Un resultado analogo al de Beale-Kato-Majda 

El lector atento se habra percatado de que la ecuacion de Burgers (15) tiene una 
especie de Teorema de Beale-Kato-Majda. Recordemos que en una dimension espacial 
lo mas parecido a u es la primera derivada espacial d x u. Entonces tenemos el siguiente 
resultado: 

Proposicion 6. Sea T > 0. Entonces, si 

fT 



/ | \d x u\ {s)ds < oo, 
Jo 



existe u solution cldsica de (15) al menos hasta tiempo T. 

Demostracion. En la prueba de la proposicion de las estimaciones a priori (Proposicion 
2.1) logramos la cota 



De esta expresion obtenemos 

— \\u\\ H 3 < c\\o x u\\l°°\\u\\ H 3, 
dt 

y aplicando la desigualdad de Gronwall (es decir, integramos la EDO) conseguimos el 
resultado. □ 



5. Buscando singular idades 

Consider amos la ecuacion (15) con 7 = 0. Esta ecuacion es una ecuacion de trans- 
pose unidimensional no-lineal, por lo tanto es susceptible de aplicarle el metodo de las 
caracteristicas. Esa tecnica aplicada a esta ecuacion se puede encontrar en casi todos los 
manuales de ecuaciones diferenciales. Nosotros utilizaremos una tecnica similar a la que 
utilizamos para probar el Principio del Maximo en la seccion 2. Queremos encontrar una 
singularidad para u. Por la Proposicion 5 tenemos que u no explota, por lo tanto si hay una 
singularidad debe estar en alguna de las derivadas espaciales de u. Como estamos en el 
caso no viscoso solamente tenemos una derivada espacial, asf que d x u es nuestra candidata 
a cantidad que explota. 

Asf sea u £ // 3 (IR) una solution clasica de (15). Por la inmersion de Sobolev se tiene que 
d x u E C 1 (M). Ahora la idea es aplicar el Teorema de Rademacher a d x u. De manera analoga 
a la de la demostracion de la Proposicion 5 se demuestra que si m x (t) = m\n x ^d x u = 
d x u(xt,t) se tiene que 

m' x (t) = d t d x u(x t ,t) = -(d x u(x t ,t)) 2 - u(x t , t)dlu(x t , t). 
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Como xt es el punto de maximo de d x u tenemos que d x u = 0, por lo que la ecuacion para 
la evolucion del maximo es 



m' = —mi. 



Si ahora resolvemos la EDO obtenemos que 

m x {t) = (23) 
Hemos obtenido el siguiente resultado 

Proposicion 7 (Singularidad). Sea f £ H ?J tal que mxiL x ^d x f(x) < 0. Entonces u 
solution cldsica de (15) con 7 = desarrolla una singularidad en tiempo finito. 



5.1. Interpretando la singularidad 

El tipo de singularidad que se produce se conoce como 'choque ' 3 y es una discontinuidad 
en u, o lo que es lo mismo d x u(x,T*) = — oo. Si pensamos en u como la altura de la 
superficie de una ola unidimensional lo que ocurre es que las particulas del fluido que 
integran la ola se mueven mas rapido cuanto mas altas esten, por lo que las mas altas van 
a acercarse a las mas bajas. Entonces cuando las alcancen se produce este choque. Esta 
discontinuidad previene que nuestra interfase u deje de ser un grafo, porque si pensamos 
en la fisica de la ola tenemos que las particulas mas rapidas sobrepasan a las mas lentas 
(ver Figura 3). 



T=0.15 T = 0.3 




Figura 3: Simulacion de una ola. 

Hemos mencionado que la ecuacion de Burgers podia ser entendida como una ola, lo 
que no es completamente descabellado. 

Si y = u(x, t) es la superficie del agua (entonces la ola viene dada por la curva 
(x,u(x,t))), y v = (v\,V2) es el campo de velocidades del agua, entonces la ecuacion 
de la interfase entre el agua y el aire es 

dtu = —d x uvi(x, u(x, t),t) + V2(x, u(x, t),t) = (—d x u, 1) • v. 
3 Del ingles 'shock'. 
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Figura 4: Una ola y sus parametros caracteristicos 

Ahora observamos que en las olas la parte se mueve a mayor velocidad que la 

parte mas baja. Por lo tanto parece natural (al menos como primera aproximacion) hacer 
la hipotesis v\{x, u(x, t),t) = cu(x, t) (la constante c podemos tomarla igual a 1 si hacemos 
un cambio de variables). Asf si despreciamos la velocidad vertical, V2(x, u(x, t),t) « 0, tene- 
mos la ecuacion de Burgers no viscosa como modelo del perfil de una ola. Podemos razonar 
tambien que la curva cambia segun sea su curvatura, de manera que V2(x, u(x, t),t) « k[u] 
con k[u] la curvatura de la curva u. Ademas si suponemos curvas suaves de amplitud pe- 
quena podemos linealizar la curvatura obteniendo el operador jd x u(x, t). Asf, con hipotesis 
razonables al menos en primera aproximacion, obtenemos la ecuacion de Burgers viscosa 
como modelo simplificado del perfil de una ola. 

Una ola unidimensional con unas ciertas caracteristicas (longitud de onda, amplitud...) 
puede modelizarse en un mejor nivel de aproximacion con la ecuacion KdV siguiente: 

dtu + d x uu = d x u. (24) 

Veamos de manera muy resumida como puede derivarse formalmente esta ecuacion. Para 
describir una ola consideramos que el agua bajo la superfie tiene un flujo irrotacional, i.e. 
v = V</> para cierta funcion escalar <f>. Si suponemos validas las ecuaciones de Euler para 
el agua bajo la superficie tenemos que (ft sigue la ley 

M+\W4>\ 2 +p + gy = 0. 

Ademas, por la incompresibilidad se tiene A</> = 0. 

La coordenada y se distingue de la x en que actua la gravedad, por lo tanto parece 
natural hacer un desarrollo de 4> en potencias de y, (j> = Y^=aV n( t ) n{x,t). Si consideramos 
olas pequehas en amplitud con respecto a la longitud de onda, entonces tenemos que 
despreciar los terminos de orden grande en y en nuestra expresion para 4>. 

Si ademas suponemos que d x (j)o w u podemos concluir la ecuacion (24). Esta hipotesis 
se motiva por los desarrollos en serie de potencias anteriores. 

Observamos que la hipotesis para obtener la ecuacion de Korteveg-de Vries es menos 
restrictiva que para obtener la ecuacion de Burgers, pues exclusivamente suponemos que 
d x <Po{x,t) = u(x,t) no que d x (j>(x , f (x , t) , t) = v±(x, f(x, t), t) = u(x,t). Consideramos por 
lo tanto discrepancias en los ordenes mayores. 

6. Metodo numerico 

En la siguiente section nos proponemos dar una aproximacion numerica al problema 
(14) propuesto en T. 
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T = 0.2 



T = 0.4 




Figura 5: Solucion para u(0, x) = — sin(x) sin viscosidad. 

6.1. Discretizacion espacial 

La forma natural de abordar el problema es usando la transformada de Fourier discreta 
(DFT) (ver [6] y [23]). Sea T^- una discretizacion de T, con un numero par N de puntos, 
dada por, 

Xj = ^(2j-N), j = 0,...,N-l 

La idea es sustituir el valor de u(-,x) = u{x) en los nodos, u(xj), por su interpolante 
trigonometrico, 

N/2-1 

I N u(x) = ^2 u k e ikx 

k=-N/2 

lNu(xj) = u(xj), j = 0, ... ,N - 1 

Es decir, dado un vector con los valores nodales de u, U = [uq ■ ■ ■ M/v-l] > l a DFT nos 
proporciona U = [u_at/2 ■ ■ ■ ^jv/2-i]- Estos coeficientes vienen dados por, 



j=0 

y puede demostrarse que estan relacionados con los coeficientes u k de la serie de Fourier 
de u(x) mediante, 

oo 

u k = u k + ^2 u k+Nm , k = —N/2, ... , N/2 — 1 

r«7^0,m=— oo 

Sabiendo como aproximar u(x), es importante saber como aproximar d x u(x). Para ello 
derivamos el interpolante trigonometrico anterior obteniendo, 



N/2-l 

£ 

k =-N/2 



8 x Inu(x) = ^2 iku k e lkx 
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Nuestra aproximacion de d x u{x) en xi vendra dada por lo que nosotros llamaremos, 

N/2-l 

(V N u)i= iku k e- ikx ', l = 0,...,N-l 

k=-N/2 

De la misma forma, parece apropiado aproximar el operador A a = {—d x ) a / 2 mediante, 

N/2-1 

Yl \ik\ a Uke~ lkxi , l = 0,...,N-l. 

k=-N/2 



T=0.75 T=1.5 




Figura 6: Solution para u(0, x) = — sin(x) con viscosidad. 

6.2. Discretizacion temporal 

Para la aproximacion de la derivada temporal en la ecuacion utilizaremos el metodo 
de Runge-Kutta explicito de orden 4, tambien denominado RK4. 

Sea T* = {to, . . . ,tjy[} una discretizacion del dominio temporal [0, T]. Asf si U s = 
[u(xi, t s ) . . . u(xn, t s )] es la aproximacion de la solution a nuestro problema en tiempo 
t = t s se tiene que U s verifica la ecuacion 

jU s = F(U s ,t s ), 

con F = —ud x u — r yA a u. 
Entonces, 

U s+ i = U S + \bt{K x + 2K 2 + 2K 3 + K 4 ), 

donde, 

K x = F(U s ,t s ), 
K 2 = F{U s + \MK 1: t s + \M) : 
K 3 = F{U S + ±AtK 2 ,t s + | At), 
K 4 = F(U S + AtK 3 ,t 8 + At). 
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